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Introduction. External photon beam modulation using compensators in order to achieve a desired dose distribution 
when brachytherapy treatment is followed by external beam radiation is a well-established technique. A compensa- 
tor modulates the central part of the beam, and the dose beneath the thickest part of the compensator is delivered 
mostly by scattered, low energy photons. A two-dimensional detector with a good spatial resolution is needed for 
the verification of those beams. In this work, the influence of different types of detectors on the measured modulated 
dose distributions was examined. 

Materials and methods. Dosimetric verification was performed using X-Omat V, Eastman Kodak radiographic films 
at different depths in a solid water phantom. The film measurements were compared with those made by ionization 
chambers. Photon beams were also modelled using EGSnrc Monte Carlo algorithm to explain the measured results. 
Results. Monte Carlo calculated over-response of the film under the thickest part of the compensator was over 15%, 
which was confirmed by measurements. The magnitude of over-response could be associated with changes in the 
spectra of photon energy in the beam. 

Conclusions. The radiographic film can be used for the dosimetry of compensated high energy photon beams, with 
limitations in volumes where photon spectra are hardly degraded. 
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Introduction 

Intracavitary application of brachytherapy sources 
followed by external beam radiation is a common 
practice in radiotherapy of carcinoma of the cer- 
vix. Since the application of brachytherapy sources 
results in characteristic dose distributions, modu- 
lated external photon beams should be added in a 
way to achieve the desired cumulative dose distri- 
bution over the target volume. Several techniques 
used in practice have been described. 14 

On the other hand, the dosimetry of modulated 
linear accelerator's photon beams is rather com- 
plex, mainly due to dose distribution in homo- 
geneity within the radiation field with large dose 
gradients. Therefore, dosimetric verification needs 



a high spatial resolution and this demand makes 
the radiographic film a dosimeter of choice. 58 
Nevertheless, it has been shown that the absorbed 
dose could be related to the energy absorption 
coefficient. 9 Furthermore, energy absorption coef- 
ficients for film emulsions and water differ signifi- 
cantly in the low energy region bellow 400 keV, 
as the data calculated according to Seltzer have 
shown. 10 Since there is an enhanced contribution of 
scattered radiation to the total dose in modulated 
photon beams, a disadvantage of the film which 
shows over-response to low energy photons may 
become important. 1113 

In this work, we compared the measured dose 
distributions of high energy photon beams ac- 
quired by different detectors. In order to discuss ex- 
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FIGURE 1. Simulated geometry of accelerator's head according to manufacturer's data for modelling compensated beam using 
BEAMnrc program package. 



perimental results, Monte Carlo simulation of par- 
ticle transport for the measured beams was done. 

Materials and methods 

Dosimetric verification of open and modu- 
lated 6 MV photon beams from Siemens Oncor 
Impression linear accelerator was performed using 
X-Omat V (Eastman Kodak) radiographic films at 
different depths in solid water (PTW Solid Water 
Phantom). We used fixed source-to-surface (SSD) 



geometry with SSD=100 cm on the phantom sur- 
face. Film dosimetry was performed using Vidar 
DosimetryPro Advantage scanner with Coherence 
Physicist (Siemens Medical Solutions) and PIPSPro 
(Standard Imaging) software packages for film 
dosimetry. The dose profiles measured by the 
film were compared with those made by ioniza- 
tion chambers (IBA Dosimetry, compact cham- 
bers CC13 and CC 01) in the water phantom (IBA 
Dosimetry, Blue Phantom). Regarding a better spa- 
tial resolution of a small volume ionisation cham- 
ber (CC01), data measured with those chambers in 



Radiol Oncol 201 1; 45(4): 310-314. 



312 



Jurkovic S et al. / Dosimetric verification of compensated beams 




FIGURE 2. Compensator's shape calculated to conform the dose distribution given by 
an external beam according to the dose distribution around brachytherapy sources. 
The insert shows thickness of the compensator in a form of level curves in mm. 
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FIGURE 3. Dose profiles measured with ionization chamber (black) and X-Omat V 
film (grey). Measurements were done with SSD=l00 cm, 20x20 cm2 field size at 10 
cm depth for open and modulated beams. Regarding the symmetry of the dose 
distributions only half-profiles are shown. 
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FIGURE 4. Mean energy distributions for open and modulated beams in air and on 
10 cm depth in water. Calculations were done with SSD=l00 cm and 20x20 cm2 
field size for open and modulated beams. Regarding the symmetry of the energy 
distributions only half-profiles are shown. 



the high gradient region of the beams were super- 
imposed on measurements with CC13 ionisation 
chamber which had a better signal-to-noise ratio. 

Photon spectra for Siemens Oncor Impression 
linear accelerator photon beams were calculated 
at the measuring planes using Monte Carlo simu- 
lation of particle transport (EGSnrc). The simula- 
tion for 6 MV photons with a field size of 2O20cm 2 
(defined at SSD=100 cm) was performed using 
OMEGA/BEAM code, developed by the National 
Research Council of Canada (NRCC). This is an 
EGSnrc user code capable of complex linear ac- 
celerator geometric coding. 14 The detailed geom- 
etry and composition of each individual device in 
the Siemens Oncor Impression linear accelerator 
were obtained from the manufacturer. Open and 
modulated beams were modelled using BEAMnrc 
software. The modelled geometry of compensated 
beam is shown in Figure 1. Calculated shape of the 
compensator is shown in Figure 2. 

The compensator shape was calculated using 
the dose distribution around the brachytherapy 
sources as a pattern according to which the open 
photon beam was modulated in order to achieve 
the desired total dose distribution. 4 

Therefore, the shape of the dose distribution is 
rather characteristic and, from the dose profiles 
point of view, three different areas can be distin- 
guished: the area under the compensator and the 
open beam area where the measurements can be 
performed with high reproducibility and the area 
near the edge of the compensator which is charac- 
terized by high dose gradients and a lower level 
of measurement reproducibility. In the last area, 
Monte Carlo calculation is especially used as a 
guideline for the interpretation of the measured 
dose distributions. On the other hand, changes in 
the energy spectrum were expected in the area un- 
der the thickest part of the compensator. 

The absorbed dose in a material depends on en- 
ergy absorption coefficients 9 and there is a large 
difference in those coefficients for film and water 
in low energy area. 10 Therefore, an over-response 
of the film under the thickest part of the compen- 
sator was expected at larger depths because the 
Compton scattered low-energy photons dominate 
there. 8 Dose calculations in a material were per- 
formed according to: 9 

D = -— = IV ^ en ^ dE = H J ^ en [1] 
P J E P P 

The BEAM code was implemented using vari- 
ance reduction techniques: photon forcing, brems- 
strahlung splitting and range rejection to speed up 
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the simulation. The lower charged particle cutoff 
energy, AE, was 0.7 MeV, and the lower photon 
cutoff energy, AP, was 0.01 MeV. The energy loss 
per transport step of the electron, ESTEPE, was 
controlled by PRESTA. 15 Scored plane was set at 
Z=100cm to collect the particles after transporta- 
tion from the accelerator, and to form the phase 
space file. Information concerning particles in the 
phase space file included the position (X, Y, Z), di- 
rection (U, V, W), energy, charge, weighting, and 
origin (LATCH). Five to ten million particles were 
collected in the scored plane. The phase space 
file served as the source for the following water 
phantom simulation using DOSXYZ, an EGSnrc 
user code for 3D absorbed dose calculation in 
Cartesian coordinates. 16 In DOSXYZ, the water 
phantom size was 40><40x40cm 3 and the phase 
space source position was on the water surface 
(Z=0). The origin was at the centre of the radiation 
field. Voxels with size of 0.5x0.5x0.5 cm 3 (XxYxZ) 
were set at the depth of the maximum dose for 
dose profile simulation. 50 voxels from water sur- 
face (Z=0) with size 2.0x2.0x0.2 cm 3 and 20 voxels 
with size 2.0x2.0x0.4 cm 3 were set along the cen- 
tral axis for central percent depth dose (PDD) 
simulation. The particles in the phase space file 
were redistributed and reused to obtain better ac- 
curacy in dose calculation. 16 Physical parameters 
of original electron beam that may influence the 
dose profile and central-axis PDD curve are the 
beam energy, the beam spot size and the distance 
from the point source. 1718 These parameters were 
adjusted to allow dose profiles and percentage 
depth dose curve to match measured data. Since 
we calculated changes in beam energies, for the 
purpose of our work, the accuracy of the beam 
profiles was not essential. We decided that 3% dis- 
crepancy from measurements is acceptable in the 
high dose region and 20% in the low dose region. 
Recommended values are 2% and 20% respective- 
ly. 19 " 21 

Results 

From the analysis of measured beam profiles, we 
observed significant discrepancies between meas- 
urements with the radiographic film and ioniza- 
tion chambers when measuring beam profiles of 
modulated beams on larger depths in water. The 
discrepancies were pronounced under the thickest 
part of the compensator (Figure 3). 

Calculated mean energy distributions in open 
and modulated beams are shown in Figure 4. 
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FIGURE 5A. Photon spectral distributions for open and modu- 
lated 6MV photon beams on central axis in air. Calculations 
were done with SSD=100 cm and 20x20 cm2 field size. 
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FIGURE 5B. Photon spectral distributions in water for open and 
modulated 6MV photon beams, on 10 cm depth in water. 
Calculations were done with SSD=100 cm and 20x20 cm2 field 
size at central axis for open beam and also at central axis 
and under steep part of the compensator (6.5 cm off axis) for 
modulated beam. 



From the 'in air 7 simulation analyses, we can 
see that the compensator removed low energy 
photons from the beam, so the mean energy of the 
modulated beam is higher than the one of the open 
beam (Figure 4). This can also be seen in Figure 5 A. 
On the other hand, at larger depths in water, the 
Compton scattering low-energy photons domi- 
nate, especially under the thickest part of the com- 
pensator, so the Figure 4 shows the decrease of the 
mean energy of the modulated beam there. 

Calculated photon spectral distributions for 
open and modulated beams are shown in Figures 
5A and 5B, respectively. 

Taking into account the dependence of mass 
absorption coefficients on photon energy for used 
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dosimeters and calculated energy distribution of 
photons in small volumes, we can estimate chang- 
es in the film response. Regarding data shown in 
Figures 4, 5A and 5B, it follows that the largest 
differences could be expected under the thickest 
part of the compensator at larger depths in water 
because of the largest energy degradation. Dose 
calculations for the film and water were done ac- 
cording to Equation [1]. Calculated over-response 
of the film in this region was over 15% which was 
confirmed by measurements (Figure 3). 

Discussion 

In this paper, we showed that magnitude of over- 
response of the radiographic film of modulated 
high energy photon beams could be associated 
with the changes in the spectra of photon energy 
in the beam. Since the largest spectral change was 
under the thickest part of the compensator, there 
was the largest difference between film and ionisa- 
tion chamber measurements. 

Regarding a high dose gradient beneath the 
steep part of the compensator, it was not pos- 
sible to measure doses in this area accurately. 
Nevertheless, spectra in this area resemble open 
rather than modulated beam beneath the thickest 
part of the compensator. In this way, the over-re- 
sponse of the film under the steep part of the com- 
pensator would be small. 

The radiographic film is often used for verifica- 
tions of modulated photon beams. 5 ' 7 ' 8 22 Despite of 
described limitations, it can be used either on build 
up depth for the evaluation of compensators shape 
or for measuring dose distributions of modulated 
high energy photon beams in phantoms. Special 
attention should be paid to the interpretation of 
measured values in volumes where photon spectra 
are hardly degraded. 



5. Bucciolini M, Band F, Casati M. Verification of IMRT fields by film dosimetry. 
Med Phys 2004; 31: 161-8. 

6. Bailey DW, Kumaraswamy L, Podgorsak MB. A fully electronic intensity- 
modulated radiation therapy quality assurance (IMRT QA) process imple- 
mented in a network comprised of independent treatment planning, record 
and verify, and delivery systems. Radiol Oncol 2010; 44: 124-30. 

7. Salz H, Wiezorek T, Scheithauer M, Schwedas M, Beck J, Wendt TG. IMRT 
with Compensators for Head-and-Neck Cancers. Strahlenter Onkol 2005; 
181: 665-72. 

8. Srivastava RP, Wagter De C. The value of EDR2 film dosimetry in compen- 
sator-based intensity modulated radiation therapy. Phys Med Biol 2007; 
52: 449-57. 

9. Brahme A, Andreo P. Dosimetry and quality specification of high energy 
photon beams. Acta Rad One 1986; 25: 213-23. 

10. Seltzer SM. Calculation of photon mass energy-transfer and mass energy- 
absorption coefficients. Radiat Res 1993; 136: 147-70. 

11. Yeo IJ, Wang CC, Burch SE. A filtration method for improving film dosimetry 
in photon radiation therapy. Med Phys 1997; 24: 1943-53. 

12. Wang X, Spirou S, Lo Sasso T, Stein J, Chui C, Mohan R. Dosimetric verifica- 
tion of intensity-modulated fields. Med Physl996; 23: 317-27. 

13. Ju SG, Ahn YC, Huh SJ, Yeo IJ. Film dosimetry for intensity modulated radia- 
tion therapy:Dosimetric evaluation. Med Phys 2002; 29: 351-55. 

14. Rogers DWO, Faddegon BA, Ding GX, Ma CM, We J, Mackie TR. BEAM: a 
Monte Carlo code to simulate radiotherapy treatment units. Med Phys 
1995; 22: 503-24. 

15. Bielayew AF, Rogers DWO. PRESTA: the parameter reduced electron-step 
transport algorithm for electron Monte Carlo transport. Nucl Instrum 
Methods B 1987; 18: 165-81. 

16. Ma CM, Reckwerdt P, Holmes M, Rogers DWO, Geiser B. DOSXYZ users 
manual. National Research, Council of Canada Report PIRS-0509B; 1995. 

17. Lin SY, Chu TC, Lin JP. Monte Carlo simulation of a clinical linear accelerator. 
App Rad Isotopes 2001; 55: 759-65. 

18. Faj D, Vrtar M, Krajina Z, Jurkovic S, Margaretic D. Model of total skin 
electron treatment using the 'six-dual-field' technique. Coll Antropol 2003; 
27: 713-21. 

19. Verhaegen F, Seuntjens J. Monte Carlo modeling of external radiotherapy 
photon beams. Phys Med Biol 2003; 48: R107-64. 

20. Rivest DRC, Riauka TA, Murtha AD, Fallone BG. Dosimetric implications 
of two registration based patient positioning methods in prostate image 
guided radiation therapy (IGRT). Radiol Oncol 2009; 43: 203-12. 

21. Mesbahia A, Reillyb Al, Thwaites Dl. Development and commissioning of a 
Monte Carlo photon beam model for Varian Clinac 2100EX linear accelera- 
tor. Applied Radiation 2006; 64: 656-62. 

22. Ting JY, Davis L. Dose verification for patients undergoing IMRT. Med Dosim 
2001; 26: 205-13. 



References 

1. Han IH, Malviya V, Chuba P, Orton CG, Devi S, Deppe G, et al. Multifractioned 
HDR brachytherapy with concomitant daily teletherapy for cervical cancer 
Gynecol Oncol 1996; 63: 71-7. 

2. Popple RA, Kim RY, Pareek P, Duan J, Shen S, Brezovich IA. Custom step 
wedge blocking using dynamic multileaf collimation for parametrial pelvic 
boost irradiation following brachytherapy carcinoma of the cervix. Med 
Phys 2003; 30: 2699-702. 

3. Jurkovic S, Zauhar G, Bistrovic M, Faj D, Kaliman Z, Smilovic Radojcic D. An 
alternative approach to compensators design for photon beams used in 
radiotherapy. Nucl Instrum Methods A 2007; 580: 530-3. 

4. Jurkovic S, Zauhar G, Faj D, Smilovic Radojcic D, Svabic M, Radiation therapy 
photon beams dose conformation according to dose distribution around 
intracavitary applied brachytherapy sources. Med Dos 2010; 35: 49-52. 



Radiol Oncol 201 1; 45(4): 310-314. 



